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We have identified a missing long-distance (one-photon exchange) contribution to the B± —> 
K ± /n£ + 1~ decays. Although it does not help to explain the anomaly measured by the LHCb 
Coll, in the ratio Rk = B(B ± —> K ± p + p~) /B{B ± —> K ± e + e~) ~ 0.75 (2.6 a deviation with 
respect to Rk = 1 + D(10~ 4 ) in the Standard Model), it provides a sizable contribution to 
the branching ratio of the B ± —» tPP t~ decays. The new decay mechanism gives rise to a 
measurable CP asymmetry, which is of order 1% in both channels. These predictions can be 
tested in forthcoming LHCb measurements. All the details can be seen in 1 . 


1 Introduction 

In the search of a more fundamental description of interactions in nature, one has to look for 
effects that can not be described by the current theory of particle physics, namely the Standard 
Model (SM). A way to do this is studying the more suppressed processes in the SM, where a 
deviation from a precise SM prediction is expected to be more significant than in processes that 
are not suppressed. In order to claim that an effect comes from physics Beyond the SM (BSM) 
in a certain process, all the SM contribution have to be known and understood at the precision 
required by experiments. In B± —> K±£ + £- decays, LHCtr has measured the Rk ratio in the 
[1, 6] GeV 2 region of squared lepton pair invariant mass q 2 , which is Rk = 0.745^qq 74 ± 0.036, 
while the current SM prediction in the same energy range i^ R ^ 1 = 1 + (S.OIq^) x 10 -4 , 
which takes into account only Short Distance (SD) contributions to the process. This LHCb 
measurement, together with a previous prediction 4 of lepton universality violation in another 
process due to purely kinematic effects, prompted us to study Long Distance (LD) contributions 
to the Branching Fraction of B± —> P±£ + l~ decays, with P = n, K. 


2 Short Distance contribution 

In order to compute the SD amplitude to the B~ —> K~/ir~l + £~ decay, we follow previous 
computations of such amplitude 5 using QCD Factorization (QCDF), where the effective weak 




Figure 1 - New long distance contribution to the process, where the square is for structure dependent one-photon 
exchange. Other structure dependent contributions vanish due to gauge invariance. 




Figure 2 - Our form factor compared to data and BaBar fitted form factor for P = n (left) and P = K (right). 
The form factors depend on the dilepton pair invariant mass. 


hamiltonian 

Heff = ~^ V ^ V ^'E C i^° q M (!) 

^ i 

is used to compute this amplitude, where q' = c,t and q = d,s. The main contributions to this 
process come from the operators with the vector and axial lepton current Og = (q'y^bi)(£ r y^£) 
and Of 0 = ((Thi^L) (^7^75-0, nevertheless at higher order in QCDF all the other operators will 
contribute, but less significantly than the former operators. 

The leading contribution to the amplitude computed from this hamiltonian is the following 


M[B- P~£ + £~] = ^(W + f aH 75*) > (2) 

V 27T 

where is the soft form factor associated to the B —>• P transition, and the Fy and Fj\ are form 
factors that depend mainly on SD effects, which becomes evident when their dependence on the 
so called Wilson Coefficients (C\, ...,C\q) is noticed. The q 2 dependence of the soft form factor 
can be obtained, for example, using Light Cone Sum Rules (LCSR) (for more details se^), 
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1 - q 2 / (5.32 GeV) 2 “ 1 - g 2 /(6.18 GeV) 2 + ^ ) 
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(3b) 


where it is fitted to dipole expressions, which are the main contributions to this form factor. 
The Pp(q 2 ) are polynomials on q 2 , whose coefficients correspond to the Gegenbauer moments. 
The values of Cg and Cio are taken from the NNLL calculation". 

































Figure 3 - New LD invariant mass spectrum from kinematic threshold for n (left) and K (right). In both plots, 
the spectrum of e + e _ invariant mass overlaps with the p + spectrum when q 2 > 0.3 GeV 2 ; therefore, there is 
a negligible contribution to Rk- 


3 New Long Distance contribution 


We focus on the same energy range LHCb uses, which excludes completely any charm quark 
contamination. The contributing Feynman diagrams are those of Fig.l; other contributions to 
this process from purely LD effects vanish due to electric charge conservation. The left-hand- 
side diagram on Fig.l is suppressed by a factor m 2 P /m 2 B with respect to the right-hand-side 
diagram, so that we can neglect its contribution. The LD contributions for the process are 
obtained using Resonance Chiral Theory (R\T) 9 . 

The amplitude that we obtain from the leading contribution is 

2 2 

M ld = V2G F V ub V: q f B f P e -i ™ B 2 [F P (q 2 ) - (4) 

H q z m B — m p 

which we see that has the same structure as the Fy part of the SD amplitude, so that it can be 
added as a correction to this form factor in the following way 


T? e ff _ zt - , K P m2 B F P(<1 2 ) ~ 1 

v ’ 9 2 M? 2 ) 


(5) 


y \f* r r 

where the kp = — 8-7T 2 , 22 , is a dimensionless constant 0(1CP 2 ) times a CKM sup- 

^ VtbV tq (mg—rrip) v K 

pression factor which depends on the final P state. In Wolfenstein’s parametrization, for P = n, 

k is 0( A°), while for P = I\ is 0( A 2 ). 


Our ignorance of the underlying dynamics is encoded in the Fp(q 2 ) form factors, which 
can be seen in Fig.2. These are almost saturated by the lowest-lying light-flavor resonances (p, 
(jj and </>), following different parametrizations foi 91 P = n and 12 P = K. We have also used 
phenomenological models by BaBai 13 as a test of our error. Their form factors are fitted using 
also heavier resonances. The almost perfect agreement on the peaks of the lightest resonances 
ensures a prediction with small error because the contribution of heavier resonances is negligible. 

The invariant mass spectrum for both channels are shown in Fig.3, which shows this spectrum 
from kinematic threshold. The difference between l = e, p becomes important at q 2 < 0.3 GeV 2 . 
We confirm that the LHCb range is free of hadronic pollution for P = K, as shown in Table 1; 
but for P = 7T there is a significant pollution in the [1,8] GeV 2 range. Comparing the LHCb 
measurement 14 to the SD contribution of the branching fraction in the whole kinematic range 1 ' 5 
we find a better agreement by adding our LD contribution, obtaining a value of BR LD+SD = 





















P — IT 

P = 7T 

P = K 


0.05 < q 2 < 8 GeV 2 

1 < q 2 < 8 GeV 2 

1 < q 2 < 6 GeV 2 

LD 

Interf 

SD 

(9.16 ±0.15) • 10 -9 
(-2.62 ±0.13) • 10~ 9 
(9.83l^) • 10” 9 

(5.47 ±0.05) • 10 -1 ° 
(-2+ 2 ) • 10” 10 
(8-71^) • 10- 9 

(1.70 ±0.21) • 10 -9 
(-6 ±2) • 10” 11 
(1.901 6 ; 6 ?) • IQ” 7 


Table 1: LD, SD and their interference contributions to the branching ratio for both channels 



q 2 > 1 GeV 2 

q 2 > 4m 2 

±3 S3 
II II 

(2.5 ±1.5)% 
-(1.3 ±0.5)% 

(14 ± 2)% 
-(0.5 ±0.5)% 


Table 2: CP Asymmetry for the different energy ranges for 7r and K. 


(2.61q 3 ) x 10 8 . The LD contribution induces lepton universality deviations of 0(10 5 ) in Rp. 
The different weak and strong phases of SD and LD contributions generate a CP asymmetry 1,16 
Acp = r(^+^p+;+;-j + r(g-^p-;+;~) i corresponding numerical results are shown in Table 2. 

4 Conclusions 

Our analysis shows that BSM studies should be restricted to the [1,8] GeV 2 range for P = tt; 
the effect of this new LD contribution could be measured in LHCb in the next run. It is also 
important to understand the current LHCb measurement of the branching fraction. In the case 
of P = K, there is not an important contribution for the branching fraction; this is because 
in the interference the peak of the (j) resonance does not surpass the CKM suppression factor, 
contrary to the pure LD contribution. The Acp we predict is an important effect that must be 
taken into account for BSM searches in both channels through this observable. 
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